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We investigate the phenomenon of the diffraction of charged particles by thin material targets 
using the method of the de Broglie-Bohm quantum trajectories. The particle wave function can 
be modeled as a sum of two terms "0 = V'mgoing + ''P outgoing- A thin separator exists between 
the domains of prevalence of the ingoing and outgoing wavefunction terms. The structure of 
the quantum-mechanical currents in the neighborhood of the separator implies the formation of 
an array of quantum vortices. The flow structure around each vortex displays a characteristic 
pattern called 'nodal point - X point complex'. The X point gives rise to stable and unstable 
manifolds. We find the scaling laws characterizing a nodal point-X point complex by a local 
perturbation theory around the nodal point. We then analyze the dynamical role of vortices in 
the emergence of the diffraction pattern. In particular, we demonstrate the abrupt deflections, 
along the direction of the unstable manifold, of the quantum trajectories approaching an X-point 
along its stable manifold. Theoretical results are compared to numerical simulations of quantum 
trajectories. We finally calculate the times of flight of particles following quantum trajectories 
from the source to detectors placed at various scattering angles 0, and thereby propose an 
experimental test of the de Broglie - Bohm formalism. 
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1. Introduction 

The aim of the present paper is t o implement the method of quantum trajectories, known as de Broglie 
- Bohm theory De Broglie , 19251 : Bohm , 19521 ] in the phenomenon of the diffraction of charged particles 
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by thin material targets. In the de Broghe - Bohm theory we consider trajectories tracing the quantum 
currents. A trajectory is determined by the initial particle's position and by the pilot wave equation of 
motion 



dr h ,Vij{r,t) 

— = — Iml — -7 ^ 

at m V'(r,i) 



a) 



where tp is the wavefunction ('pilot wave'), m is the particle mass and h is Planck's constant. This equation 
implies Newton's equation of motion in a potential 

U{r,t) = V{r,t) + Q{r,t) (2) 

where Q{r,t) is the 'quantum potential', an extra term caused by the wavefunction -0 : 



Q{r,t) 



2m 



(3) 



The probability density of an ensemble of particles guided by the same ip - field is p{r,t) = \ip{r,t)\'^. 
The equations of motion ([T]) imply the continuity equation for p. Furthermore, in the one-particle case, 
the Bohmian trajectories are equivalent to the stream lines defined by the quantum probability current 
j = (?^-/2mi')(7^*V■0— -^y-^ *). Thus, the Boh mian method yields practically equivalent results to Madelung's 
quantum hydrodynamics Madelun^ . 1926l |. 

Bohm's trajectories have been proved useful in i) constructing efficient nuinerical schemes for the in- 
tegra tion of Schrodinger's equation via swarms of trajectories (e.g. Wyattl . l2005l : ISanz et all I2nn2l : lOriolsl . 
and ii) visualizing processes for which no other intuitive picture could be obtaine d. Examples of the 



latter are the quantum tunneling eff'ect Lopreorj . 1999(|. the (particle) two-slit experim ent Philippidis et al 

ballistic transp ort through 'quantum wires' [Beenakker k, van Houtenl . Il99ll |. mole cular dynamics 
ICindens perger . 2003l |. dynamics in nonlinear sys tems with classica l focal points or caustics [Z hao &: Makri . 
I2OO3.] and rotational or atom-surface scattering |Sanz et all l2004l | . A final app lication regards the relation 
of chaotic quantum motions to the dynarnical o r igin of 'quantum relaxation' Vallentini &: and Westman 



200.^ : lEfthvmiopoulos fc and Contopoulosl . 120061 : iBennettJ . I2OIOI : IColin fc Struvvel . l201fll |. namely the justi 
fication of the approach (as time increases) of a quantum system to Born's rule p = IV'P, even if the initial 
conditions were allowed to deviate from this rule {piniUai 7^ \'4'initiai\'^) ■ 

In the case of the diffraction of charged particles, we find some new phenomena due to the structure 
of quantum currents and to their effects on quantum trajectories. In particular, we find that the quantum 
fiow is always characterized by the formation of an array of quantum vortices delineated along a locus called 
separator, i.e. a sharp boundary between the domains of prevalence of the ingoing and outgoing fiow. The 
local current structure in a vortex fori ns a pattern called 'nodal point - X-point complex'. This ha s been 
previously investigated theo retically in Efthymiopoulos et al. , 200?! : Contopoulos Efthymiopoulos . 2008 : 
Efthymiopoulos et al , 20091 ]. Here we adapt this analysis in the particular problem of particle diffraction, 
in order to determine the number and location of all critical points of the quantum fiow, as well as to 
determine the size of the nodal point - X-point complexes. The main result of this analysis is that the 
deflection of quantum trajectories is due to their approaches close to X-points. In fact, the trajectories 
follow the directions of the asymptotic manifolds of the X-points, and this fact alone suffices to explain the 
emergence of the whole diffraction pattern. 

Passing to a global description of the quantum trajectories, we identify a number of important dif- 
ferences between Bohmian and classical trajectories, which are summarized with the help of a schematic 
representation (Figure [I]). Four main differences are: 



a) In the case of quantum trajectories, the scattering angle increases when the impact parameter increases 
(Fig. la) (this is explained in section 4), while the opposite is true in the classical case. 

b) The defiection of Bohmian trajectories takes place for initial conditions at an 0{D) distance from the 
principal axis of the ingoing beam (z-axis), where D is the transverse quantum coherence length. A typical 
value for D in electron beams is D > 10~^m. On the contrary, for classical Rutherford scattering the 
impact parameter b is many orders of magnitude smaller ib ~ 10~^^m ). 



Quantum vortices and trajectories in particle diffraction 3 




z z 

Fig. 1. (a) Schematic representation of the quantum trajectories in the scattering or difFraction problem. Deflection takes 
place for quantum trajectories with initial conditions of an 0{D) distance away from the principal beam axis (z-axis). The 
thick dashed curves illustrate the separator between ingoing and outgoing (from the target) flow. The concentration of the 
quantum trajectories to various Bragg angles is shown in Figure [5] below, (b) Same as in (a) but in the approximation of 
classical Rutherford scattering for repelling ( — ), or attracting (■ • ■) forces. 



c) Bohmian trajectories have the same form in both cases of repelhng and attractive forces (Figure [T^), 
a behavior contrasted to the form of trajectories in the approximation of classical Rutherford scattering 

(FigHb). 

d) The most important difference regards the times of flight of particles to detectors placed at different 
scattering angles 6. The times of flight are shorter for larger impact parameters in the quantum case, 
independently of the sign of the interacting charges (while the times of flight depend on this sign in the 
classical case and they are, in general, much shorter than in the quantum case). This effect suggests an 
interesting experimental test that will be described in section 5. 



2. Wave function 

We consider a cylindrical beam of particles of mass m and charge Zie incident on a material target 
centered at the origin of the coordinate system of reference. Cylindrical coordinates are denoted by z 
(horizontal), R (transverse) and by the azimuth (f. Use is also made of spherical coordinates r = (z^+i?^)-*^/^, 
6 = tan~^{R/ z), and ip. 

The q uantum-me c hanical descriptio n of the diffraction process has been discussed extensively (see, for 
example, [Peng et al . 2004 : Pengj . 2005 ] for reviews and further references). Here, we adopt a simplified 



model which is sufficient for all practical purposes in the analysis below. For the interaction of the incident 
charged particles with any individual atom in the target, we adopt a screened Coulomb potential 

C/(r-rj) = J-M£!^^Ph^:;i£l^ (4) 

where Z is the nuclear charge, is the position of the j — th atom in the target and ro is a constant 
representing the screening range, whose value is taken of the order of the atomic size. The total potential 
felt by one incident particle is the sum of the individual potentials: 

N 

Vir) = Y,U{r-r,) (5) 

i=i 

where N is the total number of atoms effectively participating in the process of diffraction. 

In order to describe diffraction quantum-mechanically, we specify a model for the wavefunction -0 of 
diffracted particles (considering only elastic scattering). To this end, we first find the eigenfunctions (p of 
the time-independent Schrodinger equation: 

- —V^(l) + V{r)(l) = Ecj) (6) 
2m 
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where (j) is an eigenfunction corresponding to the energy value > of an incident particle, and V{r) in ^ 
is a time-averaged form of ([5]). We solve Eq.(l6]) following Born's approximation method (for \V\ « E)and 
expand 4> as 

<P = (j)0 + <j)l+<J)2 + ... (7) 

where (po = 0(1), (pi = OiV/E), (j)2 = 0{V'^ / E"^) etc. All the essential phenomena appear already when 
only the two first terms of the expansion — 00 + 01 are considered. We then find 



2 



l + O 



k <rj > 



(8) 



where i) the label k refers to any vector in Fourier space with modulus equal to k = {2mE)^/'^ /h, ii) 
kn = fcn with n = r/r, and iii) < r| >= (1/A^) X]j=i ^j- The corrections due to the 0{k < r| > /r) and 
1/rg terms will be omitted since they do not affect essentially the analysis. 

The time-dependent eigenfunctions are now given by 'i(^\i{r,t) = exp(— i^/c^f/2m)0k(r). The electron 
wavefunction can be modeled as a superposition of eigenfunctions 

V^(r, t) = [ d^k c(k)V;k(r, t) (9) 



i^ingoingir, t = 0) = ^ ^ X exp ( -777:^ + i{koz) ) , (10) 



(27r)3/2 

where c(k) are Fourier coefficients. Because of the collimation, the ingoing term of the wavefunction at 
t = can be modeled as a traveling plane wave along the z-direction times a Gaussian in the transverse 
direction with dispersion ~ D, i.e. 

1 / 

where, assuming that the particles are monoenergetic, the constant k^ yields the average momentum 
Po = ^^0) oi" kinetic energy Eq = h?kQ/2m along the z-direction. The value of the transverse quantum 
coherence length D, which turns to be a crucial parameter in the analysis below, depends on the details 
of the electron emission and collimation process. 

The Fourier coefficients c(k) for which ipingoing has the form (jlOp are given by: 

C(k) = jdv ^.ngoing{v,t = 0)j^ = ^^^^^^ CXp (^"^^ j (H) 

where a± = D~^. Substituting pT]) into ([9]), and using ([8]), we can evaluate the form of the wavefunction 
at all times t. After some algebra we find 

lp{r, t) = 1pingoing{r, t) + 1poutgoing{r , t) (12) 

where 

lPingoing{r,t) = ^— —^6 ^D^+^^t/^) ' ' > 13 

1 D Z,Ze^ ra Se//(fco;^,^) e-(^o-^^g*/^-) ( al\ 

and Seff{ko;9,ip), hereafter called the effective Fraunhofer function, denotes the following sum over all N 
atomic positions Vj = (xj,yj,Zj) in the target: 

N 



x] + y] 

exp ' 



2{D'^+iht/m) 
X exp ^ko{2zj sin^(6'/2) — xj sin 6 cos (p — yj sin 9 sin (p)^ 



(15) 
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The effective Fraunhofer function Seff{kQ;9,ip) accounts for the usual diffraction effects. In numerical 
calculations we consider, for simplicity, the case where the target has a polycrystalline structure. Then, 
the diffraction pattern becomes axisymmetric and the dependence of S^fj on 99 disappears. We then derive 
a fitting model for Seff (see Appendix A) which substitutes the sum (fT5]) in all subsequent calculations. 
This model reads: 



a 



q=0 

where 9q are Bragg angles defined by 



sm4(e/2)a2 2 sin [kpd smjO g) {6 - 6 g)/ 2] 

/ ^ ^coherent^ ^ i ■ I a \i a a 



(16) 



sin2(0g/2) = g = 1,2, ...q^ax (17) 

In Eq. (|16p . d is the target thickness, a is the distance between nearest atoms (assuming for simplicity a 
cubic unit cell, see appendix A). The first term in the r.h.s. of Eq. (ll6p describes a coherent contribution to 
the outgoing electron wave exhibiting sharp peaks at all Bragg angles 6g. The second term describes diffuse 
scattering due e.g. to thermal fluctuations or recoil effects in the target. These are modeled by the constant 
(Tq (dimension of length) which measures the amplitude of random motions of the atoms (see Appendix A) . 
Accordingly, the exponential factors in Eq. (jl6p are Debye- Waller factors. Finally, the constants C coherent 
and C diffuse are fitting constants whose values are fixed by a numerical simulation, while 5 is an arbitrary 
phase (see Appendix A). 

In all subsequent calculations, we use the set of Eqs. (|13|14|16p which completely specify the wave- 
function. In numerical examples we adopt the following set of parameter values: Zi = — 1, m = rn-e, 
= 8.877x lO^nm"^ (corresponding to electrons with energy E = 30KeV, or wavelength Aq = 7x 10^'^nm), 
D = lOOOnm (corresponding to a transverse quantum coherence length l//m), Z = 79 (gold), d = 420nm, 
a = 0.257nm, aa = 0.0086nm, Ccoherent = 0.060, Cdiffuse = 0.077 (found by numerical fitting), qmax = 40, 
5 = TT. These values are relevant to the case of electron diffraction, while they conform with a number of 
limitations imposed by our approximations to the wavefunction model. Some main limitations are: 

a) The energy E (determined by the particles' mass m and wavenumber ko) must be well above the 
atomic energy levels in the crystal (which are of the order of IKeV) and below the nuclear energies (which 
are of the order of IMeV). In our case the energy E = 30KeV is within these limits. Thus only the elastic 
term in the outgoing wavefunction needs to be considered. 

b) The use of Born's approximation is valid only sufficiently far from all atoms so that |y| << ii^. In 
the case of the potential of Eqs.(j3]) and ([5]), this is true at all points r of space with distance greater than 
a few times the size of the atoms in the target. All calculated Bohmian trajectories in the sequel lie in 
regions exceeding by far such distances. 

c) The inequality > (Ht/m) sets an upper limit in the time during which the ingoing wavefunction 
remains coherent in the transverse direction. Setting D = Iftm as above, we find t < 10~^sec, so that 
the total distance traveled by individual electrons should be less than Im, which is a quite long, hence a 
feasible value. 

Finally, we note that the representation of the beam as a plane wave in the z-direction can also be 
improved by introducing a finite 'longitudinal coherence length' I, i.e. a wavepacket approach in the z- 
direction as well. In this case, the fltting model of Appendix A is no longer va l id and a separate analysis 



must be made covering the cases D < I or D > I (work in progress [Delis et aLl . l201l[ |). 



3. Quantum vortices 

The form of the quantum currents in the model defined by Eqs. (|13|14|16p can be understood qualitatively 
by the following remarks: 
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3.1. Separator 

The ingoing and outgoing waves become equal in size if \ipingoing\ = \'4'outgoing\- Since ht/m « for all 
times of interest, one has from Eqs. (fT3l) and (fT4|) 



^^"Pr^J" ^ ^= V 4vreo 2fcgn2 sin2(^/2) ) ' ^'^^ 

The roots of Eq. ljlSp yield a separator line on the meridian plane (R, z = R/ tan 6) which acts as a delimiter 
between the domains of prevalence of the axial ingoing from the radial outgoing flow. Figure [5] corresponds 
to a numerical calculation, where the arrows indicate the local direction of the quantum flow at every point 
of the configuration space. Clearly, the transition from axial ingoing to radial outgoing flow takes place 
essentially through the separator lines (bold lines). 




-7500 -5000 -2500 2500 5000 7500 10000 
z[nm] 



Fig. 2. Structure of the quantum currents in the model of Eas. (|13ll4ll6p and parameters as specified in the text. 

Ray-like structures in the right part of Fig. ([2]) correspond to the form taken by the separator locally 
near every Bragg angle. Important details of the separator structure are not discernible at this resolution, 
appearing only when we zoom very close to one Bragg angle. From Eq. ijlSp we can see that roots exist only 
if G{9) satisfies G{9) < C{D), where C{D) is the maximum value of Re~^^ ^"^^^ , equal to D/y^. Then the 
two roots Ri{9) < D and i?2(^) > D define a lower and un upper branch of the separator. Only the upper 
branch i?2(^) is shown in Figure [21 since for most angles 9 the lower branch corresponds to values of Ri{9) 
smaller than R2{9) by many orders of magnitude. 

The above picture changes very close to one Bragg angle 9q, where Ri{9) becomes also important. This 
is because the local peaks of \Seff{6)\ imply also local peaks of the function G{9) in Eq. (fT8]) . According to 
the local value of the function G{9 = 9q) we distinguish the following two cases: 

Case I: G{9q) > G{D). In this case there is an interval of values 9a < 9 < 9'^ containing the Bragg angle 9q 
such that Eq. dlSp has no roots within it (Figure [3^). Then, the inner separator Ri{9) of FigJHb joins the 
outer separator R2{9) at the angles 9a and 9a and takes the form shown in FiglSb (in coordinates R vs. 
9). The gap between the left and right domains of prevalence of the ingoing flow corresponds to a narrow 
angular strip (hereafter called a 'channel') along which the flow is radial. 

Case II: G{9q) < G{D). In this case there are roots of both Ri{9) and R2{9) for all angles 9 surrounding 
and including 9q (FigiSt). However, Ri{9) has a local maximum and R2{9) a local minimum at 9q as well 
as on all other peaks caused by the side lobes of the Fraunhofer function. The separator then develops 
oscillations as shown in FiglSji. Furthermore, in this case there may be one or more pairs of angles {9, 9') 
where the inner, and outer, separator reaches the z-axis, and infinity, respectively. Two such pairs are 
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Fig. 3. (a) The function G{d) in a small neighborhood of the fourth Bragg angle Oq = 64 = 0.473811rad. The horizontal 
dashed line corresponds to the level value C{D) — De~^'^ — 606.53nm. The inner and outer separator curves are joined at the 
angles da = 0.4737825rad and O'a = 0.4738395rad, and a channel of radial flow is formed, of width AS ^ e'a-Oa = 5.7xl0"^rad. 
The function G{9) becomes zero at the two closest zeros of the effective Fraunhofer function near the Bragg peak at 9q = 64, 
namely at the angles Of^ = 0.473773 and 6^ — 0.47385. At these angles Ri becomes zero while R2 tends to infinity, (b) The 
form of the separator curves R{6) in the same range of angles 6 as in (a). (c,d) Same as in (a,b) but for a neighborhood of the 
Bragg angle 9q — 6^ — 0.676655. 



visible in FiglSjl. 



3.2. Nodal point - X-point complexes 

Along the separator, a large number of quantum vortices are formed. Their location is given by all points 
where the conditions i) ip = 0, and ii)V'i/' 7^ hold. Such points are called nodal points. The condition 
tp = 0, or ipingoing = —tpoutgoing yields, besides Eq. (fTB|) . a condition for the equality of phases, which takes 
the form (ignoring again 0{ht/m) terms) 

fco^tan((9/2) = 2g7r q Z . (19) 

The pair of equations (|18|19|) specifies completely the coordinates of one nodal point(R,^). We locate many 
nodal points numerically by choosing different values of q. It follows from Eq. (|19|) that the separation 
between two nearby nodal points, defined by two values q and g + 1, is of order ir/ko ~ Aq- 

The local form of the quantum flow in a very small neighborhood around one nodal point is very 
different from the general picture of the flow as given in Figl2J If we 'freeze' the time t, the instantaneous 
pattern formed by the vector field of q uantum probability cur r ent j c orresponds to a characteristic structure 
called nodal point - X-point complex Efthvmiopoulos et al . 2007 : Contopoulos &: Efthvmiopoulod . 2008 



Efthvmiopoulos et al . 2009l |. That is, close to a nodal point we find a second critical point of the flow. 



where one has j = 0. This is called an 'X-point', since it can be shown that it is always simply unstable, i.e. 
there are two real eigenvalues of the matrix of the linearized flow around X, which are one positive and one 
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Fig. 4. (a) Instantaneous form of the quantum flow around a nodal point placed at the left separator of the channel formed 
around the Bragg angle 64 (zq = 4505.7354nm, Rq — 2310.6028nm). The flow forms a 'nodal point - X-point' complex (Nodal 
point (N), X-point (X)) and follows in general one branch of the stable manifold (S) and then one branch of the unstable 
manifold (U) of the X-point. The other branches (SS,UU) are joined to form a loop, (b) Same as in (a) but for a nodal point 
placed on the right separator of the channel of 64 {zq = 4507.1939, Rq — 2310.9549nm). (c) Schematic representation of the 
quantum flow at the crossing of the channel, along with the directions of the stable and unstable manifolds of the X-points 
formed near the separator nodal points. The dashed lines represent the separator, (d) The size of some nodal point - X-point 
complexes, quantified by the distance Rx from the nodal point to the X-point, as a function of 6. The upper and lower dashed 
lines correspond to the estimates of Eq.((25]) for the domains of Bragg angles and of diffuse scattering respectively. 



negative. Accordingly, there are two opposite branches of unstable (U) and stable (S) manifolds emanating 
from X. On the other hand, the nodal point can be an attractor, center, or repello r. This determines the 
local form of the invariant manifolds U and S. It has been established theoretically [ Efthymiopoulos et al. 



[2007] that, except for a set of very small measure, most quantum trajectories avoid the nodal point, being 
instead scattered along the asymptotic directions of the manifolds of the X-point, leading to large distances 
from the nodal point - X-point complex. F u rthermore, while, in general, the motion of nodal point - X-poi nt 
complexes introduces chaos ( [Fris E Il997l : IWisniacki &: Puiald . l2005l : lEfthvmiopoulos et al\ . l2007l . l2009t |'l. 
in the present problem this effect is negligible because i) the speed of vortices is extremely small (of order 
~ h/ {komD"^) « vq), and ii) the quantum trajectories exhibit only a finite number of encounters with 
nodal point - X-point complexes, as will be shown with numerical examples below. In conclusion, the effect 
of the nodal point - X-point complexes on the trajectories can be described as a scattering process without 
recurrences. 

Figures |3^,b show two examples of such nodal point - X-point complexes, in which the central nodal 
points are located on the left side (FigH^) and right side (FiglDa) respectively of the main channel around 
the fourth Bragg angle shown in FiglSb- The stable and unstable manifolds of the X-points (found numer- 
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ically) are shown like thick curves in the same figures. In FigH^, the left branch of the stable manifold is 
almost horizontal far from the X-point, while the right branch of the unstable manifold follows a radial 
direction which coincides with the direction defined by the local value of the scattering angle 6. The other 
two branches form loops around the nodal point. If the terms ht/m in the wavefunction are ignored, the 
two branches join each other smoothly and the nodal point is a c enter. If however, all te r ms ar e taken into 
account, the nodal point is either an attractor or a repellor (see [Efthymiopoulos et al . 2009l |). Then one 
of the two branches in-spirals towards the nodal point, while the other deviates outwards after forming a 
nearly complete loop around the nodal point. In the present case, the time dependence of the wavefunction 
introduces negligible effects and the loops formed around all nodal points can be considered as practically 
closed. Then, in FigH^ we see that, due to the form of the stable and unstable manifolds, trajectories 
approaching the X-point from the left in a nearly horizontal direction are scattered by the X-point. A 
scattered trajectory may or may not form a loop around the nodal point. In either case, the trajectories 
eventually recede from the X-point along the unstable manifold, and follow asymptotically a radial direc- 
tion (upwards and to the right in FigU^). However, these directions are swapped in FigHja. The swapping 
can be understood with the help of a schematic figure (FigUJ:). 

Details on the size of a nodal point - X-point complex are found by expanding the wavefunction ip as 
well as the Bohmian equations of motion in variables u = z — zq, v = R — Rq around a nodal point {zq, Rq). 

Namely, from Eqs. (|13ll4p . the wavefunction takes the form 



1 



D 



-ik^ht/2 m 



where 



V2ti {D"^ + iht/mf''^ 



— g 2(_D^+iftt/m) 



sin2(0/2) 



(20) 



(21) 



and P = (ZiZe^m)/(87reo/i^). The prefactor in front of in Eq. (j20p is simplified in the Bohmian equations 
of motion, which take the form 

dz _ i;'^{di)' /dz) - ^l^'id^p'Jdz) dR _ ^lj'^{d^l^' /dR) - ^p'idip'JdR) 
'dt ~ 



(22) 



i^'i^i dt V'Vl 

Following Efthymiopoulos et al . 2009l |. the main characteristics of the equations of motion are found by 
the second order development of ip' around a nodal point {zq,Rq). For PSeff > (as in the numerical 
parameters above), the condition iP'(zq,Rq) = yields, if we disregard the term iht/m, 



PS, 



2(ro + zq) 



where ro = \/ R'q + Zq and = tan "^{Rq/ zq). Inserting this expression into the second order development 

of i)', 



tp'iu, V] Zq, Ro) = (aio + ^^lo)"" + (ooi + iboi)v 

+ ^("20 + iho)^^ + -(ao2 + ibo2)v'^ + (an + ibu)uv + 



(23) 



where u = z — zq, v = R — Rq, we find expressions for the coefficients aio, 6io jOoi j^oi; ^20 ,^20; ^02; ^02) 
Oil) bu (see Appendix B). The position of the X-point in the 'adi abatic approximation' (i.e. ig noring time 
variations of tp') can be estimated (see equations (4) and (11) of (Efthvmiopoulos et al . 2009]) as: 



~ Avx + Biu\ + Civ\ + Diuxvx 
~ -Aux + B2UX + C2VX + D2UXVX 



(24) 



where 



A = aoibio — oio&oi, 
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Tj ao2&io — aio6o2 ^ 002^10 — 010602 — 2aii6oi 7-, i, i 

-Dl — , ^1 — , ^1 — 002001 — aoiOo2 

J-, 001&02 — 002^01 ^ «oi^02 — «02^oi — 2aii6io 7-, > > 

-D2 — ^ , L'2 — ^ , ^2 — aioOo2 — 002010 • 

The variational matrix at {ux,vx) yielding the eigenvalues and eigenvectors at the X-point is determined 
by the same coefficients. 

In Appendix B we estimate the size of all coefficients aij,bij near and far from Bragg angles. Using 
the set of equations we then find an estimate for the quantity Rx = (uj^ + v^)^^'^, i.e. for the size of 
a nodal point - X-point complex. The final result is 

Rx = O ( ) domain of Bragg angles, (25) 



1 



Rx = O — -7j domain of diffuse scattering . 
\Dk^J 

Figure UJi shows the dependence of Rx on 9, obtained by computing numerically a sample of nodal point 
-X-point complexes along the separator at various angles 9. The upper and lower dashed lines correspond 
to the upper and lower estimates of Eq. ([25|) respectively. We note that the transition from the validity 
of one estimate to the other occurs around the angle 9 ~ O.Srad, This is nearly the angle beyond which 
the second term in the r.h.s. of Eq. (jl6p (diffuse scattering) becomes more important than the first term 
(coherent scattering). 



4. Quantum trajectories. Emergence of the diffraction pattern 
4.1. Trajectories 

The results of the previous section will now be used in order to understand the form of Bohm's trajectories 
of diffracted particles. A numerical calculation of a swarm of such trajectories is shown in Figure [5^. The 
initial conditions were chosen as zq = — 10/um, and xq taken uniformly in the interval [1.5, 3.3]/im, for 360 
trajectories. In FigJS^ a selective sample of trajectories is plotted so as to follow, at t = the distribution 
clN/dR oc Re^^ 1^ , with D = l^m, corresponding to the choice of ipingoing as in Eq. (fT3|) . The numerical 
integration was done with adjustable time step of maximum value At = O.lnh/E (where E = 30KeV) so 
as to ensure that the segment of a trajectory covered within one time step is significantly smaller than the 
wavelength Ao- 

The main feature of FigJS^ is that the trajectories with larger initial normal distance i?o from the z-axis 
are defiected to larger angles 9. This effect is in contrast to the picture of classical Rutherford scattering 
and it can be regarded as a clear manifestation of the quantum character of diffraction for a large transverse 
quantum coherence length. Such behavior of the quantum trajectories is readily accounted for by the fact 
that the average inclination of the separator in the (z, R) plane is negative. Namely, since all trajectories 
are horizontal until they encounter the separator, trajectories of larger initial Rq encounter the separator 
at a larger angle 9 satisfying Eq. ijlSp with R = iio(Fig.la). We thus have that 9 is an increasing function 
of i?o. 

Figure [Sb shows now a zoom in the domain of forward-scattered trajectories, leading to the formation of 
a diffraction pattern. The trajectories that started closer to the z-axis follow initially the general fiow due to 
i'ingoing, but they are subject to abrupt defiections at the crossing of any Bragg angle. Such defiections are 
due to the entering of the trajectories into channels of radial flow. This effect is depicted in detail in Figure 
[6^1, showing only one trajectory that suffers five consecutive defiections by passing through consecutive 
channels of radial flow before exiting to the domain of prevalence of the radial flow leading to the flrst 
Bragg angle 9q = 9i = 0.23. A comparison of the defiections at the Bragg angles 6q = ^2 and 9q = 9i is 
shown in panels (b) and (c) of the same figure, where the coordinates are locally rotated so that the vertical 
axis always coincides with the radial direction. We clearly see the channels of radial quantum fiow formed 
around each Bragg angle. The dashed vertical lines show local segments of the separator lines which mark 
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5000 10000 3000 6000 9000 12000 

z[nm] z[rad] 



Fig. 5. (a) A swarm of quantum trajectories (initial conditions: zq = —lO'^nm, Rq in the interval 1.5 x 10'^nm< Rq < 
3.3 X lO'^ nm). The plot shows a sample of trajectories selected by their values of Rq being distributed according to oc 

27r_Re~^ which corresponds to the choice of ipingoing as in Ea. (ll3|l . (b) A zoom of (a) in the region where the trajectories 
are forced to follow a diffraction pattern by crossing the channels of consecutive Bragg angles. The sharp deflections seen 
correspond to the Bragg angles 6q, q — 1, 8 (right to left). 




z[nm] xi[nm] Xi[nm] 

Fig. 6. (a) One quantum trajectory (bold), initial conditions {zq,Rq) = ( — lOOOOnm, 1865nm), which exits from the first 
Bragg angle after a number of visible consecutive encounters with the channels of subsequent Bragg angles (solid lines), 
(b) The crossing of the channel of the second Bragg angle 8q = 62 ~ 0.33345266 as viewed in locally rotated coordinates 
xi = z cos(7r/2 — 62) — -Rsin(7r/2 ~ 62), X2 ~ z sin(7r/2 — 62) + Rcos{n /2 — ^2)- Note the different scale on the axes. The arrows 
indicate the local direction of the quantum flow, which changes abruptly in a narrow zone near xi = 0, corresponding to 6 
being exactly equal to the Bragg value. The dashed vertical lines indicate the positions of the pair of closest zeros to the peak 
of the effective Fraunhofer function near 6 = 62- (c) Same as in (b) but for the crossing of the channel of the first Bragg angle 
9q = 81 = 0.23523778. In this case the trajectory never reaches the right border of the channel. 



essentially the limits of the channels in either case. The deflection of an orbit is caused at the crossing of 
the separator, and it is due to the orbit necessarily following the flow around the X-points that exist along 
the separator. The deflection causes a trajectory to follow a path nearly parallel to the radial direction, 
albeit with a small transverse angle, as in FigE)3, due to ipingoing, which still has some (small) influence in 
that region of the channel. Then there are two possibilities: i) the trajectory traverses the whole channel 
and exits from it from the side opposite to the entry, regaining almost horizontal flow afterwards (FigE)3), 
or ii) the orbit is entrained by the channel all along its length, in which case the particle exits from the 
scattering domain to infinity along the Bragg angle associated with that particular channel (FigjGj:). 

Since ipingoing is larger closer to the z-axis, the orbits that started closer to the z-axis have larger 
probability to cross more channels and end either with horizontal motion, or radial motion along the 
direction of one of the first Bragg angles. This leads to a complete stratification of the flow as shown in 
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O.i 1 1.5 2 2.5 

eJ[ra(l] 



Fig. 7. The difference T{6) — r(150°) of the arrival times of particles to detectors placed at the angles 6 and 150° respectively 
and at the same distance from the center, for the quantum trajectories of Fig[5^. The solid curve corresponds to the analytic 
estimate of Ea.(|26p. The arrows correspond to the first four Bragg angles. 



Figini On the other hand, due to attenuation effects (section 2) the channels practically disappear beyond 
some Bragg angle and the separator takes locally the form of FiglSjl, with oscillations of smaller and 
smaller amplitude. This causes a gradually smaller concentration of the exiting trajectories around the 
Bragg angles 9q with high q , leading eventually to a diffuse form of the radial outward flow at large angles 

e. 



5. Times of flight 

An important practical utility of the quantum trajectory approach regards the possibility to unambiguously 
determine the times of flight of particles, i.e. the time it takes for a particle to travel between an emitter 
surface and a detector surface. This question is of particular interest, because it is related to a well known 
open problem of quantum theory, namely the so-called 'problem of time' fse e iMiiga &: Leavens . 2000 : 
Muga et al . 2002 | for reviews). This problem stems from a theorem of Pauli jPauli . 1926], according to 
which it is not possible to properly define a self-adjoint time operator consistent with all axioms of quantum 
mechanics. This implies that the usual (Copenhagen)) formalism based on state vectors or density matrices 
is not applicable to a quantum-theoretical calculation of probabilities related to time observables (e.g. the 
distribution of arrival times to a detector or the times of flight defined as above). In fact, in standard 
quantum mechanics (i.e. in both Schrodinger's and Heisenberg's pictures), time is considered only as a 
parameter of the quantum equations of motion. 

Among variou s prop osals in the literature aiming to remedy this gap of standard quantum theory 
r |Muga fc Leavend . \200d/\ ). the Bohmian formalism offers a straightforward solution, since the time of 
flight is well-defined along the Bohmian trajectories. This is in principle subject to experimental testing, 
and one such test will be proposed below. We note in passing two other approaches on the same s ubjec t, 
namely the 'history approach', based on Feynman paths, and the Kijowski a pproach ( [Kijowski . 19741 ]). 
based on so-called 'Bohm-Aharonov r |Hartlel . 1 19881 : lYamada fc Takagil . ^99^) operators. However these 
approaches have not been so far as clearly formulated as the Bohmian formulation. 

A calculation of the times of flight in the case of the trajectories of Figj5]can be done as follows. We 
note first that, as a simple visual inspection of of FigH] shows, starting from a fixed horizontal distance 
from the left (source), the trajectories that are deflected to large angles are shorter in length than the 
trajectories deflected to small angles, provided that the final detection point is at the same distance r from 
the center. We can quantify this difference by approximating the form of the trajectories as horizontal 
up to their deflection at the outer separator and radial afterwards. The modulus of the velocity remains 
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nearly constant v ~ hko/m in both the horizontal and radial segments of the trajectory. Taking into 
account the separator equation (Eq. (ll8p ). we estimate the length, and hence the time difference between 
two trajectories deflected at two arbitrary angles 9i and 62 by a simple geometric analysis. If the separator 
is approximated by a straight line between the two angles (e.g. in the forward and backward directions 
01 = 30°, 62 = 150°) we find: 

TiOi) - T(92) « — (^ (-^o.2 + A(gi-g2)(cosgi-l) _ Ro,2{cos92-l) \ 
hko \ sin 01 sin02 / 

where the slope A is normalized to its value for the pair (^i, 62) = (30°, 150°) namely 

^ _ Ro,2 — Ro,i 
62 — 9i 

with i?o,i) j = 1) 2 calculated by Eq. p^ with 9 substituted by 9\ and 92 respectively. 

Eq. (j26p agrees well with a numerical computation of the times of flight for various angles 9 (see FigJT]) . 
The main remark is that, since Rq (x D (cf. Eq. (|18|) ). the time difference T{9i) — T{92) for two fixed angles 
is proportional to the transverse quantum coherence length D. That is, Eq. (|26p leads to the estimate 
T{9i) — T{92) = 0{D/vo). This is a relation that can be tested experimentally. In fact, it is very different 
to what is found in the case of Rutherford scattering (compare FigsH^ andfTJa), where a straightforward 
analysis yields 



i.e. the time difference T{9i) — T{92) does not depend on the transverse quantum coherence length D 
and it has a different dependence 0{1/vq), rather than 0{D/vq), on the particles' velocities. In the case 
of electrons, the classical time difference is of order 10"^*^ sec, while, assuming D ~ 10~^m the quantum 
time difference is of order 10~^^sec. We note that the latter time-resolution scale is within present-day 
experimental possibilities. However, in order that the experiment becomes fe asible, one shou ld be able to 
combine very accurate time measurements with single-electron detectors (see [Steinberg . 2008 ] for a review 



of recent experimental techniques on quantum-mechanical time measurements). Furthermore, one should 
also be able to have a control signal for the emission times of electrons. This is a largely une xplored subject 



'one possibility is probably offered by the so-called LASER induced field emission; see [Barwick et al. 



20071 ]). At any rate, in view of theoretical results like the above, we think it is safe to anticipate that the 



advent of experimental time measurement techniques in systems with a genuinely quantum behavior will 
open a new window for probing quantum mechanics at a very fundamental level. 



6. Conclusions 

We applied the method of the de Broglie - Bohm quantum trajectories in the problem of charged particle 
diffraction from thin material targets, focusing on the case of significant transverse quantum coherence of 
the particle beam, where new genuinely quantum phenomena appear. In particular: 

1) We constructed a model for the wavefunction of diffracted particles which takes into account both 
processes of coherent (giving rise to Bragg angles) and diffuse scattering. 

2) We developed a theory for the quantum-current structure near a locus called separator, i.e. the 
border between an inner tube of ingoing flow surrounding the beam's axis of symmetry and an outer 
domain, where the radial outward flow prevails. Analytical expressions are found for the separator. The 
separator forms thin channels of radial flow very close to every Bragg angle. Such channels are responsible 
for the concentration of quantum trajectories to particular directions of exit from the ingoing flux tube. 

3) The deflection of quantum trajectories is due to their interaction with an array of quantum vortices 
formed around a large number of nodal points located on the separator. The quantum flow near every 
nodal point takes the form of a 'nodal point - X-point complex'. The size of quantum vortices is estimated 
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analytically close to and far from Bragg angles, the estimate being in close agreement with numerical 
results. 

4) The emergence of a diffraction pattern is explained in terms of numerically calculated quantum 
trajectories. In particular, we demonstrate the sharp deflections of the trajectories as they approach one 
or more X-points along the latters' stable manifolds, and recede from these points along their unstable 
manifolds. The radically non-classical character of the quantum trajectories is demonstrated. We find that 
trajectories with larger initial distance Rq from the central axis of the ingoing flux tube (i.e. larger impact 
parameter) are deflected to larger angles 9. This is contrary to the classical Rutherford scattering, where 
the trajectories with larger impact parameter are deflected to smaller angles 9. 

5) The times of flight T{6) of particles to detectors placed at constant distances and various angles 
from the target are calculated by an analytical approximation and compared to numerical results. It is 
demonstrated that the time difference T{Oi) — T{92) follows the scaling T{9i) — T(02) = 0{D/vq), where D 
is the transverse quantum coherence length and vq the average particles' velocity. This scaling is different 
than in classical Rutherford scattering where T{9i) — T{92) = 0(|2'Zl|e^/(27^eo^T^fo)) and it leads to the 
possibility of an experimental test of the Bohmian formalism in a sector of quantum theory where the 
Copenhagen approach offers no standard recipe, i.e. the sector of time observables. 
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Appendix A Effective Praunhofer function 

We consider the simplest example of a cubic lattice arrangement of the atoms in the target. The atomic 
positions Yj are given by: 

Tj = {nx,ny,nz)a + Aauj{t), 

(A.l) 

{nx,ny,n^) £ (-^,^) x (-^,^) x (-^,^) 

where a is the lattice constant (equal to the length of one side of the primitive cell), Aa is the amplitude 
of random oscillations (due to thermal or recoil motions; Aa is taken equal to a small fraction of a) and 
Uj = {uj^x^Uj^y^Uj^z) are random variables with a uniform distribution in the intervals [—0.5,0.5]. The 
number of atoms A^^^ in the z-direction is Nz = d/a, where d is the target thickness. On the other hand, 
the value of A^j_ can be chosen as A^^ = D/a, since by such a choice the value of the Gaussian weight in 
Eq. (jl5p can be approximated by ~ 1 for all \nx\ < N±/2 and \ny\ < N±/2, and ~ for \nx\ > N±/2 or 
\ny\ > N±/2. Given these approximations, we find a model for the effective Praunhofer function (Eq. (jl5p ) 
accounting for diffraction effects when /cQfl ~ 1 (if koa >> 1 the diffraction effects disappear). In this case 
we have 

Afx/2 

S,ffiko;r)^ J2 Yl e2^^o('^-'*+^'^"^)^^°'(^/2) 

n^=-N^/2na:=-N^/2 ny=-N^/2 

^^—i[{nxa+AaUx)si-adcosLp—(nya+AaUy)si'[\ds\riip\ 2^ 

Equation ()A.2p is still not convenient for a practical use in numerical calculations, since the triple sum in 
the r.h.s. has a prohibitive cost to compute. However, a drastic simplification takes place by considering a 
material target having a polycrystalline structure. In this case, we can effectively proceed by randomizing 
the value of 99 in the terms of (IA.2P (which is mathematically equivalent to considering random rotations 
of small crystallites, on planes normal to the beam). In this case, the dependence of the sum in Eq. ()A.2p 
on Ux, Uy effectively disappears, since the corresponding terms in the exponential argument are effectively 
dominated by the random variations due to (p. Furthermore, since Uz is also random, the distribution of 



Quantum vortices and trajectories in particle diffraction 15 

values of S^jf for different samples of values {ux,Uy,Uz) becomes practically equivalent to the distribution 
of values of the quantity 

Seff {ko;r)^ e^'^" (n.a+Aau. ) sin^ (6/2) 

n,=-N,/2 

^ ^ ^ ^ ^ ^—i[{nxa+AaUx) sin 9 cos (p—{nya+Aauy) sin 8 sin {fi] (A3) 

nx=-Nj_/2ny=-Nj_/2 

i.e. where the triple sum is decomposed to a 1 x 2 sum. Denoting by Sxy an rms value of the double sum in 
()A.3p . we find that Sxy is of order A''_l D/a. We then substitute Sxy = {D /a)C coherent-, where C coherent 
is a fitting constant, in the place of the double sum. Equation (jA.Sp takes the form 

N,/2 

Seff (9) ~ (yD/a)e^^ ^ ^ Ccoherent^^P 

{i2ko{nza + AaUnJsm'^{9/2)) (A.4) 

n^=-Nj2 

where 5 is a random phase. The value of the fitting constant Ccokerent is determined numerically by running 
a realization of the sums appearing in ()A.4p in the computer, and comparing the computed quantities with 
the results with the full sum ()A.2p for various values of A^_l and Nz- To arrive at a final fitting model, we 
note first that if the noise is taken equal to zero (u„^ = 0), the form of S^ff near an angle 9q is: 

S,,eff{0) ^ {D/a)e CcoHerent k^asH9,){e - 9,) ' ^^'^^ 

Taking now int o account t he noise, the attenuation of the maxima can be estimated by a Debye- Waller 
factor (see e.g. (Pend . [2005l | 1 . changing (jA.SP into 

Sq,eff{9) ^ 

{D a)e^'Ccoherent exp --Akl sm\9 2)al ^ ■ (n \(, A.6 

V 2 J koasm[9g)(9 - 9g) 

Equation (|A.6p yields the form of the Fraunhofer function locally, very close to a Bragg angle. The total 
'coherent' contribution to the Fraunhofer function is a sum of terms like ()A.6p over all Bragg angles: 

^coherenti.^^ ~ 

fr^i ^ is''^ ( l.,2 . 2\ 2 sin [A;oiV,a sin(ej (0 - ^5)/2] 

{D/a)e'' Ccoherent exp --Akl sm\9 2)al ^, ■ \^'^ ^.7 

^ V 2 J koasm{9q){9 - 9q) 

where Qmax is the total number of Bragg angles. To this we add a 'diffuse' term accounting for random 

1/2 

phasor sums away from all Bragg angles. We have S diffuse ~ N_i_Nz , weighted by the complement of the 
Debye- Waller factor 1 — e~2^^o^^'^ (^/2)o-„_ xhus, a final model for the Fraunhofer function reads: 



Seff{0) 



sin4(e/2)a2 2 sin [/cp jV^g sin(gq) (g - gg)/2] 
4;;^°" koasm{9q){9-9q) 



q=0 



(A.8) 



where C diffuse is also a fitting constant specified numerically. Substituting Nz = d/a we obtain Eq. (fT6p . 
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Appendix B Size of quantum vortices 

Starting from Eq.([2 

Bragg angles as follows: For an angle close to a Bragg angle, using the definition (1A.5|) . one has for the 
derivatives of Seff{0) the following estimates: 



we obtain an estimate of the size of the coefficients aij, bij near and far from 



dS, 



SeffdOo 



d^S, 



eff 



Seffdel 



Niky^ikody . 



Since d/D « 1, from the above expressions we find the dominant terms of all the coefficients aij, bij. 
These are: 



OlO = e 
610 = e 

620 = e-^o/2B 



- sin(/co-zo)fco(l — cos 6*0) + 0{kQd/D) 
cos(A;ozo)^o(l — cos 6*0) + 0{kod/D) 
sin(A;o^o)^o sin^o + 0{kod/D) 
- cos{kozo)ko sin 9q + O{kod/D) 



(B.l) 



ao2 = e 



002 = e 



-Rl/2D'^ 



an = e 
611 = e-^l/^^' 



- cos(A;oZo)f sin2 Qq + Oik^/D^) 

- sin(A;oZo)f sin^ 9o + Oik^/D^) 
cos(A;oZo)f sin2 Oq + Oik^/D^) 
sin(A;oZo)f sin^ Oq + Oik^/D^) 

cos(A:oZo)f sin 2^0 + 0(^0^^7-0^) 
sin(A:oZo)f sin20o + 0(fegdVZ?2) 



(B.2) 



The quantities ro,2;o,i?o are all of order 0(D). The equalities 020 = — ao2) and 620 = —bn?. (to leading or- 
der) are due to the fulfilment of the continuity equation by the wavefunction Tp (see [Efthymiopoulos et al. 



However, since the outgoing term was evaluated only up to 0(ct^//cq) = 0(1/D^A;q), the above 
equalities are also violated at the same order, resulting in a small relative error (of order 10~^). Substi- 
tuting these expressions into Eq. (j24p . we remark that the products of the leading terms of the coefficients 
filO) ^10) 0^01; ^01 cancel exactly in the coefficient A (while they do not cancel in all other contributions). As 
a result, near a Bragg angle one has 

A = o{^Y B^ = 0{kl), a = 0{kl), Di = 0{kl), i = l,2 



whence, in view of (|24 



ux = 



Dko 



vx =0 



Dko 



(B.3) 



Substituting the above estimates in the expression Rx = (uj^ + find the first of Eqs.(l25|). 

On the other hand, far from Bragg angles all the terms with derivatives dSf>f f /dOo, d^Sf>f j /d9Q in the 
equations specifying the coefficients aij, bij become negligible. Then, we have precisely the same leading 
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terms in all the coefficients aij,bij as before. In the subsequent order, however, the terms are 0{1/D) for 
the coefficients aio,6io, aoi,&oi) and 0{kQ/D) for the coefficients 0205^207 oii,^!!, 0025^02- Then 



A = 0(^), Bi = Oikl), a = 0{kl), D, = 0{kl), i = l,2 



Thus 



which, upon substitution to Rx = (uj^ + v^)^^"^ leads to the second of the estimates ([25|) . 
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